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Abstract 

We check numerically if Soffer's inequality for quark distributions is preserved 
by next-to-leading order QCD evolution. Assuming that the inequality is sat- 
urated at a low hadronic scale we estimate the maximal transverse double spin 
asymmetry for Drell-Yan muon pair production to next-to-leading order accuracy. 
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1 Introduction 



The transversity distribution 5q(x, Q 2 ) is the only completely unknown twist- 2 
parton distribution function of the nucleon. In a transversely polarized nucleon 
it counts the number of quarks with spin parallel to the nucleon spin minus the 
number of quarks with antialigned spin [lj. In field theory the transversity dis- 
tribution is defined by the expectation value of a chiral-odd operator between 
nucleon states which is the reason why it is not experimentally accessible in 
inclusive deep inelastic lepton-nucleon scattering (DIS) |2], |3|. The most promis- 
ing hard process allowed by this chirality selection rule seems to be Drell-Yan 
dimuon production, and exactly this reaction will be utilized for attempting a 
first measurement of 5q(x,Q 2 ) at RHIC [Q. What actually will be measured is 
not the transversity distribution itself, but the transverse double spin asymmetry 
Att — d5a/da where the polarized and unpolarized hadronic cross sections are 
defined as 

d5a = - (d<r TT - da u ) , da = - (da n + da n ) . (1) 

In perturbative QCD (pQCD) Att can be expressed in terms of unpolarized 
parton distributions, the yet unknown transversity distributions and the relevant 
partonic cross sections. Although the latter have been known to next-to-leading 
order (NLO) accuracy in the strong coupling for several years by now || [| [|, 
consistent NLO calculations were not possible because of the lack of the two-loop 
transversity splitting functions. This situation changed only very recently || [| 
[T0|| , allowing for the first time for a consistent calculation of pQCD corrections 
to the transverse double spin asymmetry for the Drell-Yan process. 

The unpolarized, longitudinally and transversely polarized quark distributions 
(q, Aq, Sq) of the nucleon are expected to obey the rather interesting relation 

2\5q{x)\ < q{x) + Aq(x) (2) 

derived by Soffer (TTJ. It has been recently clarified that Soffer's inequation is 



preserved by leading order (LO) QCD evolution, i.e. if (0) is valid at some scale 



Qo it will also be valid at Q > Qq ||12|| . To NLO the situation is not as simple. 
The parton distributions are now subject to the choice of the factorization scheme 
which one may fix independently for q, Aq and Sq. One can therefore always find 
"sufficiently incompatible" schemes in which a violation of (fj) occurs. However, in 
H it was shown with analytical methods that the inequation for valence densities 
is preserved by NLO QCD evolution in a certain "Drell-Yan scheme" in which 
the NLO cross sections for dimuon production maintain their LO forms, and 
also in the MS scheme. An analytical check of the sea part is difficult since 
the singlet mixing between quarks and gluons has to be taken into account for 
the unpolarized and longitudinally polarized quantities on the right-hand-side 
(r.h.s.) of (0). In Section ^|of this article we shall show numerically that Soffer's 
inequation for sea quarks is also preserved under NLO evolution. 

Estimates of Att suffer of course from the fact that no experimental informa- 
tion on the transversity distribution is available at the moment. Therefore one 
has to rely on ansatze or model calculations of 5q(x, Q 2 ) at some reference scale 
Q fpT3|| . For example, a popular assumption is 5q(x,Ql) = Aq(x,Ql) which, 
however, is in general incompatible with Soffer's inequality (Q), in particular in 
a situation in which Aq(x, Ql) ~ — q(x, Ql)- Our aim in Sections |3] and ^ will 
be to estimate within LO and NLO an upper bound on the transverse double 
spin asymmetry for the Drell-Yan process. To do so, we will first of all assume 
validity of Soffer's inequality which seems reasonable and is corroborated by our 
finding of Section ^| that the NLO evolution to Q 2 > Ql preserves the inequation 
once it is satisfied at the input scale. The maximal asymmetry Att can then 
be estimated by further assuming saturation of the Soffer bound (||). The result 
obtained for Att under this assumption obviously strongly depends on the value 
chosen for Q . If Q is taken to be large, i.e. of the order of the invariant mass M 
of the lepton pair which sets the typical hard scale for the Drell-Yan process, the 
largest possible values for Att will be reached. However, for several reasons it 
does not seem convincing to assume saturation of (0) by the input distributions 
employing such a high Ql ~ M 2 : Firstly, evolving backwards to Q 2 < Ql - which 
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should be a completely legitimate procedure if Qo is not small - will under such 
circumstances immediately yield a violation of Soffer's inequality. Secondly, the 
r.h.s. of (HD will almost certainly lead to an overestimation for the Sq if satu- 
ration is assumed at a large Q . For instance, for sea quarks the first moment 
(x-integral) of the r.h.s. in (|2]) will diverge, which is not expected for the integral 
over Sq at any Q 2 . Therefore, to obtain a realistic estimate for an upper bound on 
A TT by assuming saturation of Soffer's inequation, two requirements have to be 
met: (i) The saturation should be adopted only at a rather low "hadronic" input 
scale where (ii) the integral over the r.h.s. of (0) is finite. Both demands are 
automatically fulfilled if we choose the unpolarized and longitudinally polarized 
input parton distributions of the radiative parton model analyses |L4|, [L5|, and 
setQ 

2Sq(x, Ql) = q(x, Qjj) + Aq(x, Q 2 ) , (3) 

where Q now is identified with the input scale \i ~ (9(0.6 GeV) of the radiative 
parton model |T3| and is considered the smallest scale from which perturbative 



evolution can be performed, such that no backward evolution from /i makes sense. 
While we are aware that our approach with its rather small Qo may lead to an 
underestimation of the maximally possible Att, we still believe our results to 



be built on a firm basis, given the large phenomenological success [|15|, [TIJ of the 
radiative parton model for the q, Aq. In any case, our results for Att under 
the assumption of (H) are the largest the radiative parton model can predict 
and will provide a useful target for future experiments. We emphasize that our 
NLO results presented in Section [| are the first ones to be obtained to true and 
consistent NLO accuracy. Section | will also provide a discussion of other possible 
uncertainties of our results. 

In Section [| we will present our conclusions. 



1 The possibility of choosing a different sign in front of the r.h.s. of (fj) will be discussed 
later. 
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2 Preservation of Soffer's inequation by NLO 
evolution 



Unlike the case of the more familiar unpolarized and longitudinally polarized 
densities, all transversity distributions obey simple non-singlet type evolution 
equations because there is no corresponding gluonic quantity due to angular mo- 



mentum conservation [17, Bf. Introducing 



5q ± (x, Q 2 ) = 5q(x, Q 2 ) ± 8q(x, Q 2 ) (4) 

and Mellin moments 5q±(Q 2 ) = J dxx n ~ 1 5q±(x, Q 2 ) the evolution equations are 
given by [|] 

Q 2 ^- 2 5q n ± (Q 2 )=5P^ ± (a s (Q 2 ))5ql(Q 2 ) . (5) 

The Mellin moments of the transverse splitting functions SPg q ± are taken to have 
the following perturbative expansion 

"SUM = (£) "2" + (S) 2 <± + • • • • (6) 

i.e. both are equal to LO. We use the following NLO expression for the strong 

coupling constant, 

a s (Q 2 ) _ 2 / A lnlng 2 /A 2 \ 

2tt /3 hiQ 2 /A 2 V Pi lnQ 2 /A 2 J ' U 

where A is the QCD scale parameter and j3 = 11 — 2rif/3, (3i = 102 — 38rij/3 

with rif being the number of active flavors. The solution of (f|) is then simply 

given by [§] 



5q n ± {Q 2 



1 + a s (Ql)-am fo^_Aj P ff 



7T/? V 2/3 

(0),n , 
Vo 

&S«3) ■ (8) 



a ,( Q2 )V 2S "«' /ft I „„ >2 



Needless to say that the LO evolutions are entailed in the above equations when 
we put the NLO quantities, £P^_j_ n , Pi, to zero. 

Eq. (H) can be very conveniently employed for a numerical calculation of the 
NLO evolution of the transversity distributions. As discussed in the introduction, 



we will assume saturation of Soffer's inequality at the input scale, see Eq. (|3]). 
Our choice for the r.h.s. of (|3|) will then be the NLO MS radiative parton model 
inputs for q(x, Ql) of |jn| and for the longitudinally polarized Aq(x, Qq) of the 
"standard" scenario of [l6j at0 Ql = [i\ LO = 0.34 GeV 2 . For simplicity we 
will slightly deviate from the actual q(x, Ql) of [O in so far as we will neglect 
the breaking of SU(2) in the input sea quark distributions originally present 
in this set. This seems reasonable as SU(2)-symmetry was also assumed for the 
Aq(x, Ql) of 0], which in that case was due to the fact that in the longitudinally 



polarized case there are no data yet that could discriminate between Am and Ad. 
We therefore prefer to assume 5u(x, Ql) = 5d(x, Ql) also for the transversity 
input. We will examine the possible effects of SU(2)-breaking later. The moments 
of the resulting input distributions 5q(x, Ql) are easily taken, and the 5q±(Ql) 
are then evolved to higher scales Q 2 > Ql with the help of (||). A standard 
inverse Mellin transformation finally gives the desired transversity distribution 
in x-space. In order to perform this inverse Mellin transformation, Eq. (§) has 
to be analytically continued to complex n [|TJ[]. The evolutions of the q(x, Ql) 
(neglecting the SU(2)-breaking) and the Aq(x, Ql), which both involve the singlet 



mixing between quarks and gluons, proceed as explained in [jTJ, [15|, [16| 



In order to numerically check the preservation of @, Fig. |T] shows the ratio 

Rq{x > Q2) = q(x } Q 2 ) + Ag(x, Q 2 ) (9) 

as a function of x for several different Q 2 values for q = u v = w_, u = — n_)/2, 
d v = gL, d = (d + — gL)/2 (cf. Eq. @). If NLO evolution preserves Soffer's 
inequality then R q (x, Q 2 ) should not become larger than 1 for any Q 2 > Ql- As 
we already know from || this is the case for the two valence distributions. Fig. [1] 
confirms that Soffer's inequality is indeed also preserved for sea distributions. 
Furthermore, in Fig. || we see that evolution leads to a strong suppression of 
R q (x,Q 2 ) at small values of x, in particular for the sea quarks. This can be 

understood by the fact that 5P qq ^±(x) has a very mild behaviour for x — » ||, 

2 Note that for the purpose of checking the preservation of Soffer's inequality by evolution 
the choice of the initial scale Ql is actually irrelevant. 



and by the well-known sharp small- a; rise of the unpolarized sea distributions in 
the denominator of R q due to Q 2 -evolution. We note that our numerical results 
for the sea quarks became somewhat unstable at large x, probably caused by the 
fact that the sea distributions are obtained here as differences of two much larger 
quantities. 

As is obvious from (0), Soffer's inequality only restricts the absolute value of 
the transversity distribution. Therefore, we are free to choose a different sign 
in front of the r.h.s. of (§) and have to check the results for the two distinct 
possibilities Sq v (x, Qq) > 0, Sq(x, Ql) > and Sq v (x, Ql) > 0, Sq(x, Q 2 ) < 0. Our 
results do not noticeably depend on the actual choice. 

As we have neglected any possible SU(2)-breaking in all the sea input dis- 
tributions q(x,Ql), Aq(x,Ql), 6q(x,Ql), any difference between the curves for 
Ru, Rd can necessarily only result from the dynamical breaking of SU(2) first 
induced by NLO evolution. The occurrence of a small breaking from this source 
is well-known from the unpolarized [18j and longitudinally polarized cases. 
For the transversity densities it is given by 

MQ 2 ) -a s (Ql)) f xp (l),n _ jj P (l),n 
Tl(3o 

' a s {Q 

.a.(Q8) 

Fig. displays the resulting effect via the ratio 



IT On V / 



xf^IT (6u v - 5d v ) n (Ql) . (10) 



, 2 _ Su(x,Q 2 ) - 8d(x,Q 2 ) 

l3W } 5u(x,Q*) + 6d(x,Q 2 ) [ ' 

for various Q 2 . One can see that - apart from the region of very large x - the 
dynamical breaking of SU(2) is rather small and could in reality well be entirely 
masked by the explicit breaking in the non-perturbative sea input. 
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3 Upper bounds on Att- Framework 



Now that we have shown that NLO evolution preserves Soffer's inequation, we 
want to utilize it to derive upper bounds on the transverse double spin asymmetry 
to be measured in polarized Drell-Yan muon pair production. For this purpose we 
choose again the maximally allowed value @ for the transversity distributions, 
which should yield the maximal double spin asymmetry. We employ the same 
unpolarized and longitudinally polarized input distributions as in the previous 
section, along with the same value for the initial scale Qq. 

The scaling variable for the Drell-Yan process is r = M 2 /S, where M is the 
invariant mass of the produced muon pair and yS is the center-of-mass energy 
of the hadronic collision. Since in unpolarized reactions only the collision axis 
is specified, the distribution of the produced muon pairs cannot depend on the 
azimuth cf). If the colliding nucleons are transversely polarized then the collision 
and spin axes specify a plane in space and consequently the polarized cross section 
will depend on 0. Instead of working with r-dependent cross sections we again 
prefer Mellin moments defined by 

t\ TT ^mL . (12) 



d(j) J Q drdcf) 
Including NLO corrections to these cross sections one obtains the generic expres 
sion [|, 



d(5)a n a 2 ei 



cm 



9S 



(<W) 



m n q {Ql) (i + {*)c? Y > n {Ql) 



+ H-(Ql)^f- (6)C? Y > n (Q* F 
where 

(5)H;(Ql) = Y,4[(^A(Ql)(m(Q 2 F ) + (A^B)] , (14) 
q 

H n g {Q 2 F) = T,<[9a(QI)^b(QI) + ^b(Qf)) + ^^B)] . (15) 
q 

The dependence on the azimuth is given by $(0) = 1 and £$(</>) = cos 20. 
Integration over <fi thus isolates the unpolarized part and $(0) is then replaced 



(13) 



by 2tt. On the other hand the integration @ (/!^/ 4 - + JsT/^ ~ IsH/i) ^ 
extracts the polarized cross section and 8<f>(<f>) can then be simply substituted 
by 4. In the following we will always assume appropriate integration over the 
azimuth. 



The unpolarized NLO MS QCD coefficients in r-space can be found, e.g., in 
|}J. Their Mellin moments are 



C° Y > n {Q 2 F ) 



C F [ ASf(n) 

c F 2 

Tr 



n(n + 1 
+ 3-4S 1 (n 



Si(n) 



+ 



rr 



{n + iy 



+ 



4vr s 



+ 



n(n + 1 

n 2 + n + 2 
n(n + l)(n + 2) 



In 



M 2 

Of 



(16) 



n 2 + n + 2 



In 



S'i(n) + 



M 5 



-ll^ + 22n 2 + 14n + 4 
n 2 (n + l) 2 (n + 2) 2 



n(n + l)(n + 2) ~ V.<2f. 
where CV = 4/3, Tr = 1/2. The polarized ones can be found in 
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(17) 



5C» Y > n (Q F ) = C F 



4Sf(n) + 12 (S 3 (n)-C(3)) 



and read 
4tt 2 " 



n(n + 1) 



C F [3-45i(n)] In 



5C^ y ' n (g|) 



(18) 
(19) 



In the above formulas we used the abbreviation Sk{n) = Y^j=ij k - Since there 



is no gluon transversity distribution for the nucleon, the gluonic part of (|I~3|) 
drops out for the polarized case. The indices A and B in ([14]) and ( 15]) take into 
account the possibility of having two different scattering hadrons, although only 
pp collisions are planned at the moment. Finally, Q F and Qr in Eqs. (O)-(O) 
are the factorization and renormalization scales, respectively, for which we will 
choose Qf — Qr = M unless stated otherwise. 

Z° production and 7Z°-interference can be easily included by the substitution 
(see also [f?]) 

2 



» M 2 (M 2 - My) 
p 2 _ <Ap V,V K - Z 

% ^ VlV ^(M 2 -M 2 z ) 2 + T 2 z M 2 z 

M 4 

+ i6 ( y 2 + a 2 ) { v; ± ^ 2 (M2 _ M|)2 + r|M| 



(20) 



S 



where 

k = ^ FMl , V f = Tf-2e f sin*e w , Af = Tj . (21) 
167ra em j j j j j 

The positive sign in front of the A 2 term is appropriate for the unpolarized cross 



section, the negative sign for the polarized one. As usual, Gp denotes the Fermi 



constant, Qw the Weinberg angle (sin 2 Qw — 0.224) and T? the third component 



of the weak isospin. 

For examining the perturbative stability of our results we will also calculate 
the cross section at LO. In this case one simply needs to set the QCD coefficients 
to zero in the above formulas and to replace the NLO parton distributions by 
ones evolved in LO. As LO input distributions for (J3|) we will use the unpolarized 
LO parametrizations of |TjJ (neglecting again the SU(2) breaking in the quark 
sea) and the polarized ones of [|16| at the LO input scale = 0.23 GeV 2 . 

4 Results 



Fig. |3| shows the transversely polarized pp cross section and the "maximal" double 
spin asymmetry App for v^S" = 40 GeV, corresponding to the proposed pT 



HERA-N fixed target experiment which would utilize the possibly forthcoming 
polarized 820 GeV proton beam at HERA on a transversely polarized target. We 
show results at both LO and NLO. For illustration we have also included the 
expected statistical errors for a measurement of App by HERA-N which can be 
estimated from 

SApp = 1 , (22) 
P B Pp\ Cat 

where Pb and Pp are the beam and target polarizations for which we will use 
Pb = Pp = 0.7. C is the anticipated integrated luminosity of C = 240 pb -1 , a the 
unpolarized cross section integrated over bins of M, and e the detection efficiency 
for which we will take for simplicity e = 100%. Note that full 4n coverage of the 
detector is assumed. Fig. ||| shows that the maximal asymmetry for HERA- N is 
actually fairly large and would be accessible in that experiment. 



In Fig. ^ we present results similar to Fig. [| but now for = 150 GeV, 
corresponding to the RHIC collider. We note that the region 9 GeV < M < 11 
GeV will presumably not be accessible experimentally since it will be dominated 
by muon pairs from bottomonium decays. Again the predicted maximal asym- 
metry is of the order of a few per cents. From the expected error bars calculated 
again for 70% beam polarization, C = 240 pb -1 and e = 100% one concludes that 
asymmetries of this size should be also well measurable at RHIC. 

Fig. |5| shows similar results for the high-energy end of RHIC, ^fS = 500 GeV, 
where the integrated luminosity is expected to be C = 800 pb -1 . It turns out that 
the asymmetries become smaller as compared to the lower energies, but thanks to 
the higher luminosity the error bars become relatively smaller as well, at least in 
the region 5 GeV < M < 25 GeV where the errors are approximately 1/10 of the 
maximal asymmetry. One can also clearly see in Fig. |5| the effect of Z-exchange 
and the Z resonance. 

We have already mentioned before that Soffer's inequation does not determine 
the sign of 5q(x, Q 2 ), so that in principle we have to check all different combina- 
tions in order to find the "true" maximal value for Att- It turns out, e.g., that 
keeping a positive sign only for 5u v (x, Q 2 ) leads to a reduction of \Att\ at small 
M but an enhancement at the experimentally not accessible region of large M. 
We have checked that for small M the asymmetry takes its largest values if all 
signs are chosen to be positive, as was done in Eq. (H) and in the above plots. 

A comparison of the LO and NLO results in Figs. §-|5| answers one key ques- 
tion concerning the transversely polarized Drell-Yan process: Our predictions for 
the maximal Att show good perturbative stability, i.e. the NLO corrections to 
the cross sections and Att are of moderate size, albeit not negligible. There 
seems to be a general tendency towards smaller corrections when the energy is 
increasing, which should be mainly due to the larger invariant masses probed and 
to a resulting smaller a s (M 2 ). 
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Let us finally address some of the the main uncertainties in our predictions 
for the maximal asymmetry Att- The first issue is the scale dependence of 
the results. This is examined in Fig. || for the case \fS = 150 GeV. We plot 
here the maximal asymmetry in LO and NLO, varying the renormalization and 
factorization scales in the range M/2 < Qf = Qr < 2M. One can see that 
already the LO asymmetry is fairly stable with respect to scale changes, which is 
in accordance with the finding of generally moderate NLO corrections. The NLO 
asymmetry even shows a significant improvement, so that Att becomes largely 
insensitive to the choice of scale. 

In order to get a rough idea about the uncertainty caused by our imperfect 
knowledge of the longitudinally polarized parton densities Aq(x, Q 2 ), Ag(x, Q 2 ), 
we have also calculated the asymmetries using the NLO "valence" scenario input 
distributions of |TB| instead of the "standard" ones in (|3p. As can be seen in 



Fig. [7| for the case \/S =150 GeV, the difference for experimentally significant 
M turns out to be quite small, with the predictions based on the "valence" 
scenario distributions having slightly smaller asymmetries. 

As we already mentioned in Sec. |2|, neither the "standard" nor the "valence" 
scenario parametrizations take into account a possible SU(2) breaking in the 
polarized sea because only neutral current polarized DIS data are available at the 
moment. This led us to neglecting also any SU(2) breaking in the transversity 
input densities for our calculations, just keeping the dynamical SU(2) breaking 
produced by NLO evolution (cf. Fig. |2|). On the other hand, it seems rather 
likely that a certain amount of SU(2) breaking - possibly much more than the 
one generated by evolution - could be realized in nature. One possible way of 
estimating the uncertainty entering our predictions for Att through this source, is 
to reintroduce the hitherto neglected amount of SU(2) breaking in the unpolarized 
input densities as fixed in the original input distributions of [|TjJ. The SU(2) 
breaking will also influence the transversity input via Eq. (|3]). The resulting 
asymmetry is also depicted in Fig. [7]. As can be seen, the effect is sizeable only 
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at rather large M. 



5 Conclusions 

We have shown numerically that Soffer's inequation is preserved by NLO QCD 
evolution provided it is satisfied by the input distributions. 

For the first time, we have presented a complete and consistent NLO calcu- 
lation of the transverse double spin asymmetry A TT for the Drell-Yan process, 
employing the NLO corrections to the hard subprocess cross sections as well as 
performing the (^-evolutions in NLO. Here we have estimated the maximally 
possible Att in the framework of the radiative parton model by assuming that 
Soffer's inequality is saturated at a low hadronic scale. For \/S = 40 GeV the 
maximal value of Att for pp collisions was found to be approximately 4% with 
an expected statistical error for HERA-N of about 1% at an invariant mass of 
M = 4 GeV. The situation for RHIC with = 150 GeV turns out to be 
rather similar. The prospects of measuring Att somewhat improve when going 
to \/S = 500 GeV where the maximal asymmetry is of the order of 1% for small 
M with an expected relative statistical error of approximately 1/10. We empha- 
size again, however, that our results only represent an upper bound on Att, so 
that the "true" asymmetry may well be much smaller and even experimentally 
not measurable. 

Comparing to corresponding LO calculations, we find that the QCD correc- 
tions turn out to be moderate but non-negligible, putting our predictions on a 
firm basis. We have also examined the main uncertainties of our predictions, 
such as the scale dependence of the asymmetry and our imperfect knowledge of 
the longitudinally polarized parton densities to be utilized for the saturation of 
Soffer's inequality at the input scale. We found that these uncertainties seem 
to have rather little impact on our results in the regions hopefully accessible in 
future experiments with transversely polarized protons. 
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Note added: After completing this work, we received the paper [2^] in which 
a mathematical proof of the preservation of Soffer's inequality under NLO Q 2 - 
evolution is given. 
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Figure 1: The ratio R g (x,Q 2 ) as defined in (Q) for q = u v ,u,d v ,d and several 
fixed values of Q 2 . 
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Figure 2: The dynamical SU(2)-breaking in the NLO transversity densities ex- 
pressed by the ratio SD(x, Q 2 ) as defined in QTTJ) for several fixed values of Q 2 . 
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Figure 3: NLO and LO maximal polarized Drell-Yan cross sections and asym- 
metries for HERA-N. The error bars have been calculated according to Eq. ( p2f) 
and are based on C = 240 pb" 1 , 70% polarisation of beam and target and 100% 
detection efficiency. 
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Figure 4: As in Fig. |, but for RHIC at = 150 GeV assuming £ = 240 pb" 1 , 
70% polarisation of each beam and 100% detection efficiency. 
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Figure 5: As in Fig. |, but for VS = 500 GeV and C = 800 pb _1 . 
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Figure 6: Scale dependence of the LO and NLO asymmetries at y/S = 150 GeV. 
The renormalization and factorization scales in (|ll^ -(|l"5"|) were chosen to be Qr = 
Q F = M/2, M, 2M. The solid line is as in Fig. § 
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Figure 7: The NLO asymmetry at y/S = 150 GeV, using the "standard" and the 
"valence" sets of |R| for the Aq(x, Qq) in (^). Also shown is the change caused 
by not neglecting the SU(2) breaking in the q(x, Qq) of (jTSfl (see text). The solid 
line is as in Fig. |j. 
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